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. Abstract 

The solutions of the reduced Maxwell-Bloch equations for an anisotropic two-level medium, which describe 
the propagation of electromagnetic pulses having a duration from a few field oscillations, are studied. An 
influence of the permanent dipole moment of the quantum transition on dynamics of the pulses and their 
■ spectrum is considered. 
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1 Introduction 

GO 

The significant attention was paid during the last years on investigation of the nonlinear coherent phenomena 
in the anisotropic media 1-10 . Distinctive feature of such the media is that the stationary states of quantum 
particles being contained in them possess no the parity. For this reason, the diagonal elements of the matrix of 
the dipole moment operator and their difference commonly referred to as permanent dipole moment (PDM) of 
t-H ' the transition are distinct from zero. 

The study of the formation and dynamics of electromagnetic pulses in the anisotropic media has great 
, practical importance. Modern technologies give an opportunity to produce the semiconductor crystals with 
properties varying over a wide range. For example, an anisotropic medium, in which the modulus of the ratio 
between PDM and corresponding off-diagonal matrix element of the dipole moment operator lies for different 
quantum transitions within the limits from 0.15 up to 7.1, is resulted in work 1 . Besides, one of the persistent 
trends in developing the laser physics is to obtain the most shorter pulses 11,12 . 

The role of PDM in the second harmonic generation in asymmetric semiconductor quantum wells was recog- 
nized and extensively studied (see, e.g., 1 and references therein). Its relative influence on the conversion efficiency 
increases for larger pump/harmonic wavelengths, especially at higher pump powers. Also, propagation in the res- 
onant optically uniaxial media of the two-component electromagnetic pulses consisting of the short-wavelength 
• -h , ordinary and long-wavelength extraordinary components was investigated in 6-9 . It was revealed that, in the 
conditions of strong interaction between the components, pulses can pass through the anisotropic medium in the 
^ , modes distinct from the self-induced transparency. 

The propagation of one-component extremely short electromagnetic pulses of femtosecond durations through 
the media having PDM was considered in works 2 " 5,10 . As opposite to the ultrashort pulses, a high-frequency 
carrier wave is absent in the Fourier spectrum of such the pulses. It was noted in 2 that nonzero PDM may 
have a strong effect on the medium response for very short pulses, while an application of the usual scheme 
of the rotating-wave approximation eliminates its influence in the longer pulses case. Complete integrability 
of the reduced Maxwell-Bloch equations for the two-level medium with PDM in the frameworks of the inverse 
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scattering transformation method 13-15 was established in paper 3 . The pulse solutions of these equations on 
constant (but not arbitrary) background were obtained. The background was chosen in such a manner that the 
problem considered was reduced to the isotropic medium case. Dynamics of these pulses at the pump presence 
was considered in 4 . The system of the full Maxwell-Bloch equations was studied in work 5 , where stationary 
extremely short pulses decaying exponentially and algebraically were found. An asymmetry caused by PDM 
on the polarity of these pulses was revealed. Also, the effects of passing the extremely short pulses with a 
duration up to several oscillations of electromagnetic field through a medium possessing PDM were investigated 
numerically 10 . It was shown that there exists solitary stable bipolar signal with a nonzero time area (nonzero 
breather). The recent results summarized above make valuable to clarify the role of PDM in the formation the 
one-component pulses of the self-induced transparency (including, the extremely short pulses and ultrashort ones 
as the limiting cases). In present report, we shall study the exact solutions of the breather type of system of the 
reduced Maxwell-Bloch equations for an anisotropic medium. 



2 The model and basic equations 

Let's consider optically uniaxial medium, whose anisotropy is created by the electric field. It splits the energy 
levels owing to the Stark effect, retaining the degeneration of electronic levels on the absolute value of projection 
M of the total angular momentum. Thus, 7r-transition (AM = 0) and doubly degenerate er-transitions (|AM| = 
1) are formed in an electronic subsystem. 

Let the electromagnetic pulse propagate in the positive direction of y axis perpendicular to the optical axis 
z of the medium. Suppose that only the extraordinary component E e of its electric field, which is parallel to z 
axis, is distinct from zero. One can show that such the pulse will interact with the 7r-transition only. The system 
of reduced Maxwell-Bloch equations describing this process in the unidirectional propagation approximation 16 
has next form: 

— =l -E e (a-a), (1) 

— =i(GJ + -E e )a + 2i-E e <T 3) (2) 
dE e n e dE e .Nduj 

— + — = -2m (cr-cr), (3) 
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where 173 = (/?22 — Pn)/2 is the population inversion of the quantum level; a = p\2', Pjk {k = 1,2) are the elements 
of the density matrix; d, D and ujq are the dipole moment, PDM and the frequency of 7r-transition, respectively; 
n e is the extraordinary refractive index; N is the concentration of the ^-transitions. In the equations given, we 
ignore the relaxation and inhomogeneous broadening of the spectral line of resonant absorption. 

System coincides with the reduced Maxwell-Bloch equations for the isotropic media 16 if D = 0. One 

can see that the extraordinary component of electric field fulfills here two functions: it causes the quantum 
transitions and, simultaneously, shifts dynamically their frequency. In the case of the two-component pulses 6-9 , 
these functions were executed by its short-wavelength ordinary and long-wavelength extraordinary components 
accordingly. In so doing, the ordinary component generated extraordinary one due to PDM of the resonant 
cr-transitions. 

It is convenient for subsequent consideration to introduce new variables 

dE e / n e \ nNd 2 
u = - — , t = lj [t y), 77 = 2 —y. 4 

Tiloq V c / n e cn 

Then, equations CEJ-© are rewritten as 

dr 



3 = iu(a — a*), (5) 



do 

— = i(l + 2ku)(T + 2iua 3 , (6) 
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du 
dr) 



where 



i(a* - a), 
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An integrability of equations JSJl 0) with arbitrary parameter k in the frameworks of the inverse scattering 
transformation method 13-15 was revealed in 3 . The corresponding overdetermined system of linear equations 
(Lax pair) has the next form 



OT 
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where ip = ip(r, rj, A) = (V>i, ^2) T is the vector solution of the Lax pair, A is so-called spectral parameter, matrices 
L(X) and A(X) are defined as follows 
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One can check by direct calculation that the compatibility condition of both the equations of Lax pair JSJ 

dL(X) dA(X) 



dr] 



[L(X),A(X)}=0, 



(9) 



is equivalent to system 



3 Exponentially and rationally decreasing breather-like pulses 



Being integrable by the inverse scattering transformation method, equations J5J— |(7J have the multisoliton solu- 
tions. The one-soliton pulse and its algebraic limit can be obtained by direct integration of these equations in 
stationary case as it was done, for example, for the system of the full Maxwell Bloch equations 5 . The formulas 
arising at it differ from ones presented in 5 by the definition of the pulse velocity only. Applying the Darboux 
transformation technique 17 to find the two-soliton breather-like pulse solution of system ©-J7J) on the zero 
background, we come to the following expression for variable u: 



d 



VTTfc 2 dr 



— arctan 



4fcarfisinh A R + Z+T^ 1 cosAr 



Z_ £1 cosh A R - AkapT- 1 sin Ai 



(10) 



where 



A R 
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VI + k 2 
CI 



(r + 16(1 + k 2 ){Z\ + AZ + 16fc 2 4)^^) 
r - 16(1 + k 2 )(Z 2 -AZ- 16fc 2 af) 



+ Ci, 



C2, 



A = Z\+ %Z\{Z - Ak 2 a 2 R ) + 16(1 - Aa 2 R ){Z 2 - 4fc 4 a|), 
fl = a R Z+/Z, T- 1 =a I Z-/Z, Z = A(a 2 R + a]) , Z ± = Z±k 2 . 



(11) 
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Real constants an, aj, c\ and C2 are the free parameters of the pulse. Constant do is an initial population of the 
quantum level (|<Jo| < 1/2). The formulas for variables 03 and a are cumbersome and omitted. 

The time area of this two-soliton solution is distinct from zero: 



4sign(fc) 
udr = ; arctan 



\Tk\ 



a? R + aj 



For this reason, we refer to it as the breather-like pulse to distinguish from the well-known breather pulse in 
isotropic media, whose time area is equal to zero. An existence of the breather type solutions with nonzero area 
was established numerically for system ©-0 in 10 . 

Quantities f2 and T, which characterize the frequency and duration of the pulse, can be regarded as its free 
parameters instead of an and a/. Then, one finds from equations 

aR = 2 [ 1 + frf 

l + R- 



2T 
where 

1 



R± = ^r±{{n 2 -k 2 )T 2 -l), 

r = {{Ft 2 - k 2 )T 2 + I) 2 + Ak 2 T 2 . 

Here the arithmetic root is taken in the definition of r, and the signs of i?+ and i?_ are chosen in such a manner 
that R+R- — FIT . Without loss of a generality, we shall suppose in the sequel that parameters an and a/ (or 
Fl and T) are positive. If FIT ^> 1, then solution 1)1 OJl is nothing but the pulse with higher- frequency filling 
(ultrashort pulse) , while at FIT < 1 variable u may no change the polarity as it happens for the extremely short 
pulses. 

Expanding the right-hand side of equation HOf) in the Taylor series at a neighborhood of point k = and 
retaining the first two terms, we come to formulas 

U = Uq + ku\, 

FIT cosh Br sin Bj + sinh B r cos Bj 
U ° = r> 2 T 2 (cosh2 J B i? , + 1) + cos IB 1 + 1 ' { ' 

4» 2 r 2 (cosh25 fl + 1)(cos2Bj + 1) 
Ul ~ ~ (n 2 T2(cosh2B R + 1) + cos 25/ + l) 2 ' ( } 

where 

r 4T((fl 2 + 1)T 2 + l)cro77 



Br = — 
Bi = fir 



T ((n + l) 2 T 2 + l)((ft-l) 2 T 2 4 
AQT 2 ((n 2 - 1)T 2 + l)aoV 



((n + i) 2 t 2 + i)((n - i) 2 t 2 + 1) ■ 

It follows from the expressions presented that the modulus of the Fourier transform of uq achieves the maxima 
on odd harmonics of the basic frequency f2, while the modulus of the Fourier transform of u\ has them on even 
harmonics, including zeroth one. Since the generation of the secondary harmonics due to PDM is nonlinear effect, 
they are more localized at the center of the pulse. In the case of breather-like pulses, an asymmetry induced 
by PDM manifests itself in that the signs of k and the zero harmonic are opposite. Similar asymmetry on the 
polarity of a signal for the extremely short pulses was revealed under analytical and numerical investigations 
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The conclusions made above are most obvious for the pulses with slowly varying envelope. Indeed, if condition 
QT ^S> 1 holds, then formulas (|12|l and (| 1 31) become simpler: 



2 sin Bi 

u = 7FZ r o ■ ( 14 ) 



T cosh Br 



2 



nr cosh b r 

Carrying out the Fourier transform with expressions l|14|l and (fL5 



, 2 cos B j 

«i = " o,^ _ • (15) 



OO 

F{v,u ,i) = J e tiyT u 0A dT, 



one obtains 

F(,, Uo ) = ^exp^T^) (^____- ______ j 

cxpH__ fl »/) ^ 2^ e X p2i(0 J -OT0 ii ) 
ui) = -TT — — + {v + 211)- 



fi 2 V sinh7rTiy / 2 y sinh7rT(^ + 2^)/2 

o ., exp2i(„__ fl -0 f ) ^ 
+ [ ) smh7rT(u-2Q,)/2j' 

where and 9j are the values of Br and 5/ at r = 0. 

It is seen that the maxima of the moduli of F(y,ua) and F(y, m) are reached at v — and at v — 0, 2S_, 
respectively. The width of the spectral lines is equal to T _1 , while the maximum values of \F(v, u±)\ are 
proportional to fl~ 2 T~ 1 . Thus, an efficiency of the secondary harmonic generation grows with reducing the 
carrier frequency of the pulse. This agrees with the results of the previous studies (see, e.g., 1 ). An influence of 
PDM on the pulses having the high-frequency filling is weak. This is due to the fact that an average value of 
the frequency detuning 2ku (see formula on the duration of such the pulse tends to zero. 

Let us return to expression (|10fl again. The time dependence of u and 173 is plotted below. Parameters 




Fig. 1. Profiles of u and 03 with k = 1, ao = —0.5, fi = 1 and T = 10. 

n = n/y/l + k 2 and f = y/l + k 2 T of the pulse are selected so that it excites the medium strongly (see Fig. lb). 
The curve of the modulus of the Fourier transform of F(i>, u) is given on Fig. 2. Note, that the zeroth harmonic 
exists in the Fourier spectrum in accordance with that the time area of the breather-like pulses differs from zero. 
For the pulses with QT w 1, this fact was established in 10 (see Fig. 3 also). Since the principal maximum on 
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Fig. 2. Modulus of Fourier transform F(u,u) of the pulse presented on Fig. 1. 




Fig. 3. Profiles of u with k = 1, (To = —0.5, O = 1, T = 1 (solid curve) and with fc = (dotted line). 



Fig. 2 is arranged to the left from O, the secondary harmonics are generated most effectively by the pulses, whose 
basic carrier frequency on an input of the anisotropic medium is lesser than the resonant frequency. The same 
effect takes place under reducing a duration of incident pulse due to broadening of its spectral width. Also, it 
is possible to show that the position of the principal maximum on v axis in the Fourier spectrum is displaced 
to a red region with growing The curve on Fig. 4 shows a difference of the spectral structure of it's with 
identical O and T in the cases of anisotropic and isotropic media. The lapse to the right of O is a consequence 
of increasing an asymmetry of the principal peak of the Fourier spectrum in the medium with PDM. 

In limiting case 2(an + £<xj) — > kenpiip (or T — * oo) formula IjlOjl yields 

17 cot 9? Ci sin 2 iy9 cos £2 + (cos 2 </? — 2) sin £2 + 2sin(^ 



y/1 + k 2 (£1 sint^ + cos £2) 2 + cotV (1 - sinv? sin £2) 2 ' 



where 

40(1 + O 2 )^?? 4Hoo?7 
Cl= ° T+ (1-02)2 +«*. C2=nr + T -^ + c 2 , 

k cos </5 

~ x/TTfc 2 ' 
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Fig. 4. Difference AF of the moduli of Fourier transform F(y, u) of the breather-like pulses with 
k — 1 and k = (isotropic case). Values of the remaining parameters are the same as for the pulse 
plotted in Fig. 1. 



Here variable u decreases rationally. An arbitrary parameter of the pulse is real constant tp determining its 
carrier frequency, which is always less resonant one. 

Existence of rationally decreasing pulses is the distinctive feature of the anisotropic media. So, algebraic 
one-component extremely short pulse was found in 5 , and two-component one-parameter pulses were constructed 
in 9 . Curves of u, 03 and \F(y, u) \ for the pulse are represented in Figs 5 and 6. Comparing curves on Figs 2 




Fig. 5. Profiles of u and 0-3 of the rationally decreasing breather-like pulse with k = 3, 00 = —0.5 
and ip = tt/8 (fi w 0.88). 

and 6, one can see that the position on v axis of maxima of the secondary harmonics is more shifted in the red 
region for the rationally decreasing pulses. Besides, the expressed asymmetry has not only the principal peak, 
but the secondary peaks of the Fourier spectrum also. 

Performing limiting procedures tp — > ±7r/2 and —c\ — > ±cosc2, we obtain from equation (|16|l : 

= 2fc 

l + k 2 + k 2 (T + 4a r]) 2 ' 

If we impose additional condition C2 — > ±7r/2 in so doing, then, instead of (|17|l . formula (|16fl turns into following 

u = -6fc— , (18) 
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\F{v,u)\ 




Fig. 6. Modulus of Fourier transform F(u, u) of the rationally decreasing breather-like pulse. Values 
of the parameters correspond to the pulse on Fig. 5. 

where 

id = k 4 (r + 4<7 J]) 4 + 6fc 2 (r + 4(7077) ((1 + k 2 )r + 4(1 - 3k 2 )a v) - 3(1 + fc 2 ) 2 , 

w 2 = k 6 (r + 4<7 Q r]) 6 + 3fc 4 (r + 4cr r/) 3 ((1 + k 2 )r + 4(1 + 9fc 2 )<r r/) + 

+ 9fc 2 (3(l + fc 2 ) 2 r 2 + 8(1 + fc 2 )(3 - k 2 )a Q Ti] + 16(3 - 2k 2 + HfcVoV) + 9(1 + fc 2 ) 3 . 

At the end of this section, we discuss briefly the spectral data 13-15 corresponding to the solutions studied. 
The discrete spectrum of breather-like pulse H1Q(I includes four simple points. They are located on the plane of 
the Lax pair © spectral parameter at 



A = ± 



A(a R + iai) 2 - k 2 



A = ± 



4(or - ia/) 5 



2{a R + iai)Vl + k 2 ' 2(a R - ia/)\/l + k 2 

The points of the discrete spectrum of rationally decreasing pulse (|16|l are doubly degenerate and located at 

ik sin if 



A = ±- 



VTTfc 2 ' 



Solution (|17|l has two simple points of the discrete spectrum at A = Hk/y/1 + k 2 . This rational pulse is a 
particular case of the one-soliton solution of system which is stationary. The points of the discrete 

spectrum of solution (|18(l are the same as for previous one, but they are doubly degenerate. 



4 Conclusion 

In the present report, we have considered the propagation through the two-level medium possessing PDM of the 
electromagnetic pulses with duration from a several oscillations of the field. The explicit breather-like solutions 
of corresponding system of the reduced Maxwell-Bloch equations are studied. Unlike to the breathers in an 
isotropic medium, these ones have the time area distinct from zero. An existence of the solutions of such a 
kind (namely, the nonzero breather) was revealed under numerical analysis in 10 . As in the case of the extremely 
short electromagnetic pulses 5,10 , it is shown here that an asymmetry on the polarity of a signal takes place 
for the breather-like pulses of arbitrary duration: the signs of the zero harmonic and PDM are opposite. It 
follows from the Fourier analysis of the solutions under discussion that an influence of PDM on an efficiency of 
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the secondary harmonics generation grows with a reduction of the basic carrier frequency of the pulse and its 
duration shortening. As a result, the pulses, whose carrier frequency on an input of the medium is less than 
resonant one, will generate the secondary harmonics most effectively. This concern, in particular, the rationally 
decreasing breather-like pulses that exist only in the case of the anisotropic media. 

5 Acknowledgment 

This work was supported by the RBRF grant 05-02-16422 a. 

References 

[1] S. Kocinac, Z. Ikonic, and V. Milanovic, "The influence of permanent dipole moments on second harmonic 
generation in asymmetric semiconductor quantum wells", Opt. Communic. 140, 89-92, 1997. 

[2] L. W. Casperson, "Few-cycle pulses in two-level medium", Phys. Rev. A57, 609-621, 1998. 

[3] M. Agrotis, N. M. Ercolani, S. A. Glasgow, and J. V. Moloney, "Complete intcgrability of the reduced 
Maxwell-Bloch equations with permanent dipole", Physica D138, 134-162, 2000. 

[4] A. A. Zabolotskii, "Pumping of the extremely short pulses in optical medium", JETP 94, 869-882, 2002. 

[5] A. I. Maimistov, and J.-G. Caputo, "Extremely short electromagnetic pulses in resonant media possessing 
the permanent dipole moment", Opt. Spectrosc. 94, 245-249, 2003. 

[6] S. V. Sazonov, "Resonant transparency effects in an anisotropic medium with a constant dipole moment", 
JETP 97, 722-737, 2003. 

[7] S. V. Sazonov, and N. V. Ustinov, "Optical transparency effects in the modes of long-short-wave length 
resonance", Proc. SPIE 5402, 72-80, 2003; "The effects of optical transparency in the modes of resonance 
of long and short waves", Izv. Ross. Akad. Nauk, Scr. Fiz. 68, 1280-1283, 2004. 

[8] S. V. Sazonov, and N. V. Ustinov, "Resonant transparency regimes under conditions of long/short- wave 
coupling", JETP 100, 256-271, 2005. 

[9] S. V. Sazonov, and N. V. Ustinov, "The optical transparency regimes in anisotropic media", Izv. Ross. 
Akad. Nauk, Ser. Fiz. 69, 1136-1138, 2005. 

[10] S. O. Elyutin, "Dynamics of an extremely short pulse in a Stark medium", JETP 101, 11-21, 2005. 

[11] D. H. Sutter, et al, "Sub-6-fs pulses from a SESAM-assisted Kerr- lens modelocked Ti: sapphire laser: At 
the frontiers of ultrashort pulse generation", Appl. Phys. B: Lasers Opt. 70, S5-S12, 2000. 

[12] V. Loyo-Maldonaldo, et al, "Generation of ultrashort electrical pulses in semiconductor waveguides", IEEE 
Photon. Technol. Lett. 15, 428-430, 2003. 

[13] S. P. Novikov, S. V. Manakov, L. P. Pitaevsky, and V. E. Zakharov, Theory of Solitons: the Inverse 
Scattering Method, Consultants Bureau, New York, 1984. 

[14] G. L. Lamb, Jr., Elements of Soliton Theory, Wiley, New York, 1980. 

[15] Solitons, Edited by R. K. Bullough, and P. J. Caudrey, Springer, Berlin, 1980. 

[16] J. C. Eilbeck, J. D. Gibbon, P. J. Caudrey, and R. K. Bullough, "Solitons in nonlinear optics I. A more 
accurate description of the 2n pulse in self-induced transparency", J. Phys. A: Math. & Gen. 6, 1337-1347, 
1973. 

[17] V. B. Matveev, and M. A. Salle, Darboux transformations and solitons, Springer- Verlag, Berlin-Heidelberg, 
1991. 



9 



